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Abstract 



By a symbolic method, we introduce multivariate Bernoulli and Euler polynomials as 
powers of polynomials whose coefficients involve multivariate Levy processes. Many properties 
of these polynomials are stated straightforwardly thanks to this representation, which could 
be easily implemented in any symbolic manipulation system. A very simple relation between 
these two families of multivariate polynomials is provided. 

keywords: multivariate moment, multivariate Bernoulli polynomial, multivariate Euler poly- 
nomial, multivariate Levy process, umbral calculus. 

1 Introduction 

Quite recently many authors have obtained a moment representation for various families of polyno- 
mials. For the multivariate Hermite polynomials H v (x), with v = (v i, . . . , vj) € Ng a multi-index, 
i.e. a vector of nonnegative integers, the moment representation 



has been given by Withers [TT] with E the expectation symbol, i the imaginary unit, Y ~ 
N(0, and £ a covariance matrix of full rank d. Making use of the Laplace distribution and 
of the Gamma distribution, Sun [TIT] gives a moment representation for Bernoulli polynomials, 
Euler polynomials and Gegenbauer polynomials in the univariate case; see also references therein. 

By using a symbolic method, known in the literature as the classical umbral calculus [7], a 
different moment representation for multivariate Hermite polynomials is provided in [I], without 
using the imaginary unit. Umbral methods are essentially based on a symbolic device consisting 
in dealing with sequences of numbers, indexed by nonnegative integers, where the subscripts are 
treated as powers. Under suitable hypothesis (see [2] for a detailed discussion), these sequences 
of numbers could be interpreted as moments of random variables (r.v.'s). 

In this paper we show how the classical umbral calculus allows us to give a moment represen- 
tation like (jl.ip also for the multivariate Bernoulli polynomials B^\x) and the Euler polynomials 
£v(x), where an additional real parameter t € R is included. While in (jl.lj) the random "part" is 
represented by the multivariate Gaussian random vector Y, in the representation here introduced, 
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the random "part" is represented by a multivariate Levy process [8]. Thanks to this representa- 
tion, we point out many similarities and a new relation between these two families of polynomials. 
Open questions are addressed at the end of the paper. 

2 Multivariate umbral calculus 

The classical umbral calculus has reached a more advanced level compared with the notions that 
we resume in this section. We only recall terminology, notation and basic definitions strictly 
necessary to deal with the topic of the paper. We skip the proofs, the reader interested in is 
referred to pjj 2] . 

A univariate umbral calculus consists of an alphabet A = {a, f3, . . .} of umbrae, and an eval- 
uation linear functional E : R[A] i— > R, defined on the polynomial ring M.[A] such that 

i) E[l] = 1; 

ii) E[a l (3 3 ■ ■ ■ ] = Ela^Elfi 3 } ■ ■ ■ (uncorrelation property) for distinct umbrae a, (3, . . . and non- 

negative integers 

A sequence a$ = l,ai,a2, ... S 1 is umbrally represented by an umbra a if i?[a n ] = a n , for all 
nonnegative integers n. The element a n is the n-th moment of the umbra a. The same sequence 
of moments could be represented by infinitely many and distinct umbrae. More precisely, the 
umbrae a and 7 are said to be similar if = £7 [7"] for all nonnegative integers n, in symbols 

a = 7. 

An umbra looks like the framework of a r.v. with no reference to any probability space. The 
way to recognize the umbra corresponding to a r.v. is to characterize the sequence of moments 
{a n }. When this sequence exists, we can compare the moment generating function (m.g.f.) of the 
r.v. with the so-called generating function (g.f.) of the umbra, that is the formal power series 

f(a,z) = l + Y,a n — (2.1) 

n>l 

If the moments of the r.v. are defined only up to some finite m, then one works with sequences 
of m elements only; see |9j for further details. 

Let us consider a <i-tuple of umbral monomials /x = (/ii, . . . , /i^) and set fx" = fi^ 1 ■ ■ ■ /j^. A 
sequence {g v } v&N d G R, with 

9v — 9vi,...,v d & n d go — 1; is umbrally represented by the d-tuple \i if 
E[/j, v ] = g v , for all v G Ng. The elements {g v } ve j^d are the multivariate moments of fj,. If {^i}f =1 
are umbral monomials with disjoint supports then g v = E[^] ■ ■ ■ E\/i v /]. The g.f. of the d-tuple 
fx is 

z) = E[e^+-+^] = l + J2 E 9v ~^\ ( 2 - 2 ) 

k>i ven^,\v\=k 

where z° = z^ 1 ■ ■ ■ z^ d , \v\ = v\ + ■ ■ ■ + Vd and v\ = v\\ ■ ■ ■ v^l. Two d-tuples /x x and /x 2 are said to 
be similar, in symbols /x 1 = // 2 , if and only if f(fi l ,z) = f(fi 2 i z )i that is ^[/^i] = Elfx"} f° r all 
v G Ng. They are said to be uncorrected if and only if Elfiifj,™] = E^WE^^] for all v,w G Nq. 

Multivariate Bernoulli umbra. Let t be the Bernoulli umbra [7], that is the umbra with 
g.f. f{i,z) = z/(e z — 1), whose n-th coefficient is the n-th Bernoulli number. 

lr The support of an umbral polynomial p £ R[_4] is the set of all umbrae in A which occur in p. 
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Definition 2.1. The multivariate Bernoulli umbra t is the d-tuple with all elements 

equal to the Bernoulli umbra l. 



From Definition 12.11 and (|2.2p . we have 



f(t,z) = E[e 



tziH Viz d -\ 



f(t,zi H V z d ) 



Definition 2.2. The multivariate Bernoulli numbers {B^ } vE ^d are the coefficients of the g.f. 
IP) , that is Bi 1] = E[l v }. 

Since E[t°] = E[l Vi l V2 ■ ■ ■ i Vd ] the following result is proved. 



Zi-\ h z d 

e 2iH Yz d _ i ' 



(2.3) 



Proposition 2.3. B^ ] = E[i\ v \) for all v £ 



o- 



(fci) ' ' ' (k d ) ^ or v = • • • ' v d)i k = (ki, ■ ■ ■ ,kd) £ Nq, and assume k < v if and 
only if fcj < f j for all j G {1, 2, . . . , d}. 



Set (2) 



Proposition 2.4. j = £ fc<1> ' for all DeN|j suc/i that \v\ > 1. 

Proof. Let w be the e?-tuple (it, . . . , u) with all elements equal to the unity umbra it, that is the 
umbra with all moments equal to 1. We have 



e(;K w =e(;H*i=e(;H^ 

fc<« v 7 fc<« v 7 fc<« v 7 



and by linearity pQ 



k<v ^ 7 



fc<1) 



L k U V ~ k 



E[{l + u) v ] = E[{l + u 



.Hi 



as £?[(t + = £?[nj=x(i + u) Vl \- Since + = E[(t + l) k ] = E[i k ] for all nonnegative 
k > 1 [7], the result follows from Proposition 12.31 for k = \v\. □ 

Multivariate Euler umbra. Let r\ be the Euler umbra, that is the umbra with g.f. f(r], z) = 
2e z /[e 2z + 1], whose n-th coefficient is the n-th Euler number. 

Definition 2.5. The multivariate Euler umbra r] is the d-tuple (?/,••• , rj), with all elements equal 
to the Euler umbra r\. 



From Definition 12.51 and (|2,2p . we have 

f( V , z) = E[e^ + - + ^} = f(n, z x + ■ ■ ■ + z d ) 



2 e ( z iH — 

e 2(z 1 + ---+z d ) _|_ i " 



(2.4) 



Definition 2.6. The multivariate Euler numbers {^v } ve ^d are the coefficients of the g.f. 
that is = E[r) v \. 

Proposition 2.7. <*i 1} = £[r/H] f or a u v e pjg. 
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Multivariate Levy processes. One feature of the classical umbral calculus is the feasibility 
to extend the alphabet A by adding new symbols [7J, the so-called auxiliary umbrae, whose 
moments depend on moments of elements in A. A very important example is the so called dot- 
product m.a of a nonnegative integer m and an umbra a. By using the exponential Bell polynomials 
[3], the moments of m.a can be expressed in terms of moments of a, since m.a represents a 
sum of m uncorrelated umbrae similar to a. So we have f(m.a,z) = [f(a,z)] m and similarly pQ 
/(m./Lt, z) = [/(/x, z)] m . Thanks to the notion of auxiliary umbrae, in this last equality the integer 
m could be replaced by a real number t € K 

f(t.iM,z) = [f(»,z)] t . (2.5) 

Indeed the multivariate moments of m./x are pQ 

E[{m.nY] = Y, w^TiM/(A)^[MA] (2-6) 
^m(A)!A! 

where the sum is over all partitions A = (A^ 1 , \ 2 2 , . . .) of the multi-index E[[i^] = g^g^ • • • , 
with <7,\i multivariate moments of /x, and (jn.)& denotes the lower factorial. Prom (|2.6p . -E'[(m./x)' u ] 
results to be a polynomial q v (m) of degree |u| in m. Then, the symbol t.fj, denotes the auxiliary 
umbra such that E[(t.fj,) v ] = q v (t), by which (|2.5p follows (see Proposition 2.2 in [1]). In particular, 
for the multivariate Bernoulli and Euler umbrae we have 

/»•*■«> =( £^ + ji )' "'■■"'-( ..^i )' < 2 - 7 > 

The auxiliary umbrae i.t and i.r/ are symbolic versions of multivariate Levy processes. Indeed, 
let X = {X t }t>o be a Levy process on M d , that is a stochastic process starting from and 
with stationary and independent (i-dimensional increments. According to the multivariate Levy- 
Khintchine formula [8], if we assume that X has a convergent m.g.f. tp x (z) in some neighborhood 
of 0, then we have 

Vx(z) = (<p Xl (z))\ (2.8) 

with X\ = (x[ , . . . , Xj ). Within the multivariate umbral calculus, if we denote by /x the d- 
tuple such that /(/x, z) = (p Xl (z), then the auxiliary umbra t.fi is the umbral counterpart of X. 
The auxiliary umbra t./x has various algebraic properties paralleling those of t.a pQ. We recall 
those we will use later on: 

t.(cn) = c(t.fi), (t + s)./x = i./x + s./x, t-(/Xi + /x 2 ) = t.fi 1 + t./x 2 (2.9) 

for c, s,i € M, with s ^ t, and /x x and fx 2 uncorrelated d-tuples of umbral monomials. If we 
replace s with —t in the second equivalence of (j2.9j) . the auxiliary umbra — i./x has the remarkable 
property — t./x + t./x = e, where e is an umbra with g.f. /(e, z) = 1. The auxiliary umbra — i./x is 
called the inverse of i./x and it is such that —i./x = i.(— l./x). Then also —i./x is a symbolic version 
of a multivariate Levy process. As example, to keep the length of the paper within bounds but 
also for the open questions addressed in the last section, we just show the probabilistic counterpart 

2 A partition of a multi-index v, in symbols A h v, is a matrix A = (Xij) of nonnegative integers and with no 
zero columns in lexicographic order such that A r i + X r 2 + • • • + Xrk = v r for r = 1, 2, . . . , d. The length /(A) of A 
is the number of columns of A. The notation A = (Ap , , • • •) means that in the matrix A there are ri columns 
equal to Ai, V2 columns equal to A2 and so on, with Ai < A2 < • • ■ . We set rn(A) = (n, V2, ■ ■ ■)■ 
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of —t.L and —t.rj. H Indeed the following propositions give the probabilistic interpretation of the 
d-dimensional random vector X\ in (|2.8p corresponding to — l.i and — l.Tj respectively. 

Proposition 2.8. The inverse —l.i of the multivariate Bernoulli umbra is the umbral counterpart 
of a d-tuple identically distributed to (U, . . . , U), where U is a uniform r.v. on the interval (0, 1). 

Proposition 2.9. The inverse —l.tj of the multivariate Euler umbra is the umbral counterpart 
of a d-tuple identically distributed to (X, . . . ,X), where X = 2Y — 1 with Y a Bernoulli r.v. of 
parameter 1/2. 

Definition 2.10. The t-th-order multivariate Bernoulli numbers {B$ } ve ^d are the multivariate 
moments of the multivariate umbra t.L, that is B§ = E[(t.L) v ]. 

Definition 12.101 generalizes the definition of the multivariate Bernoulli numbers given in [5]. 
In particular we have B^ = E[(0.i)°] = 1 and s£ 0) = E[(0.i) v ] = if [v| > 0. 

Proposition 2.11. B® = ^ k <v il) ^ ^v-k > f or al1 s,teRandve Ng. 

Proof. For s = t, the proof is straightforward. For s,t € R with s ^ t, from the second of (j2.9|) . 
we have t.L = (t - s).l + s.l, so that E[(t.L) v ] = E[{(t - s).i + s.l} v ] for all v € N$. The result 
follows from Definition EEfl since E[{(t -s).l + s.l} v ] = J2 k <v E [{(t ~ s)-i} v ' k ]E[(s.L) k }. □ 

In Proposition 12.111 set t = 0. We have the following corollary. 

Corollary 2.12. For all s 6l we have J2k<v (fe) E k v^k = ^ if v = otherwise being 0. 

Definition 2.13. The t-th-order multivariate Euler numbers {&v} V £f$d are the multivariate mo- 
ments of the multivariate umbra t.rj, that is = E[(t.r]) v ]. 

Definition 12.131 generalizes the definition of the multivariate Euler numbers given in [5] . As 
before we have = ^[(O.ry) ] = 1, and ^ = E[(0.rj) v ] = if |v| > 0. Proposition EH] and 
Corollary [2712] can be restated also for <£$ since their proofs depend only on the moment umbral 
representation and not on the properties of the involved umbrae. 

3 Multivariate Bernoulli and Euler polynomials 

In the classical umbral calculus, we can replace the field R with R[xi, . . . , Xd], where x±, . . . ,Xd are 
indeterminates [T]. Then the linear operator E is defined on the polynomial ring R[xi, . . . , Xd] [A] 
with values in R[xi, . . . , xj. The only hypothesis to be added on the linear operator E is that if 
x = (xi,...,x d ) then E[x v ft w ] = x v E[/j. w ], for all v,w G Ng. 

Definition 3.1. (Moment representation of multivariate Bernoulli polynomials) The 

multivariate Bernoulli polynomial of order v 6 Nq is B^\x) = E[(x + t.i) v ], where l is the 
multivariate Bernoulli umbra and t G R. 

Probabilistic counterparts of t.L and t.r) could be given, but the involved random variables are less known. This 
goal goes beyond the aim of the paper. 



5 



Definition 3.2. (Moment representation of multivariate Euler polynomials) The mul- 
tivariate Euler polynomial of order v E Nq is £v\x) = E {(a; + \ \t.{r) — u)\) v \ with igR,« = 
(u,...,u) a d-tuple with all elements equal to the unity umbra u and rj the multivariate Euler 
umbra. 

Definition 13.11 and 13.21 means that the multivariate Bernoulli polynomials and the multivariate 
Euler polynomials are such that 

„v—k 



k<v 



k<v 



V \ X 



k 2l fe l 



-E[{t.( V -u)} k ] 



(3.1) 



where t.t and t.{j] — u) are multivariate Levy processes. This symbolic moment representation 
of the coefficients simplifies the calculus and is computationally efficient [HE]- In particular, the 
definition of multivariate Euler polynomials is given according to the terminology first introduced 
by Norlund [6]. Indeed, from the first equivalence in (|2.9[) we have \[t.(r) — u)] = t. — u )} ■ 
Thanks to Proposition ^. 91 we have —1. [^(?7 — w)] = i [ — 1./7 + u] and this last symbol represents 
a (i-tuple identically distributed to (Y, . . . , Y), where Y is a Bernoulli r.v. with parameter 1/2. 

Proposition 3.3. B$(x) = Zk<v {Z) xV ~ kB k and 2^8® {\x + |l) = J2k<v Q* v ~ k ^k with 
1 the d-tuple with all elements equal to 1. 

Proof. The former equality follows from the first equality in (|3.1|) . by observing that E[(t.i) v ] = 
By ■ For the latter equality, observe that 



E 



Since E 



1 t 1 r , 



k<v 



E 



t 1, s 
-1 - -(t.u) 
2 2 V ' 



[x + t.T]) 



v—k 



2\v-h\ 



{|l — ^(t.ii)} fc = for all fc, except when k = which gives 1, we have 
2M£® (±x + ±l)=E [(x + t.rj) v ] = g Q x"- k E[(t.T,) k ] 



by which the result follows. 



□ 



Corollary 3.4. B$ = B^ (0) and <E$ = 2^1 (|l) , mi/t and 1 i/ie d-tuples with all 
elements equal to and 1 respectively. 

Corollary 3.5. E[B$\-t.i)\ = E[B$ (t.(-l.i))] = and E[£® [\{t.(u - l.rj)])} = 0. 

Taking into account Definitions 13. 1[ 13.21 and (|2.7p . the g.f. of the multivariate Bernoulli and 
Euler polynomials are respectively 



f(x + t.L, z) = e 



xizi-\ Yx d z d 



glj Kgrf 

e 2iH \-z d _ 1 



,f [ x + -[t.(r] -u)],z 



2tgZiziH \-XdZd 



Proposition 3.6. If x = (x\, . . . ,Xd) and y = (y±, 
then 



(e z i+-+ z d + 1)' 
, yd) are two d-tuples of indeterminates, 



B^\x + y) = £ ( V M\x)B^_ k (y) £^\x + y) = £ ( V M\x)S<>l k (v). 



k<v 



k<v 
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Proof. We replace M[xi, . . . , x^\ with W[xi, . . . , x&, Hi, ■ ■ ■ ,yd\- From the second equivalence in 
(|2.9p . we have 

B { * +S \x + y) = E{[(x + t.t) + (y + s.l)] v } = £ (YW* + t.L) k ]E[(y + a.Cf-% 

k<v ^ ' 

by which the former equality follows. The latter equality follows by the same arguments. □ 
Corollary 3.7. E k < v ®4\x)si-l<x) = E k < v ©^V)^-^) = 

Proof. In Proposition 13.61 se t s = —t. The result follows by observing that (2x) = B^\2x) = 
E[(2x) v ] = 2\ v \x v . □ 

Proposition 3.8. E$(tl - x) = (-1)1*1^ (x), and £ v t] {tl - x) = (-1)^6^ (x). 

Proof. The former equality follows by observing that Bv\tl — x) = E[(t.u — x + t.t) v ] and 

E[(t.u-x + t.t,) v ] = (-l)^E[(t.(-u) + x + L(-l)) v ] = (-IpE {[x + t.(-(t + u))] v }. 

Since E[(i+u) k ] = (-l) k E[i k ] for all nonnegative integers k [7], then E[(-(l+u)) v ] = (-1)H.E[(M- 
u)\ v \] = E[l\ v \\ = E[l v }. Tl 
follows. Similarly we have 



= E[iW] = E[i v ]. Then we have -(t + u) = i and £.(-(fc + w)) = t.t, by which the result 



(tl —x) = E 



tl-x + t.- - -1 
2 2 



(-i)Hf; 



Since /(— ^,z) = the latter result follows. □ 

As it happens for the univariate case, the multivariate Bernoulli and Euler polynomials share 
many properties. Undoubtedly, this is due to the connection between Bernoulli and Euler num- 
bers that here is emphasized by the similar multivariate moment representation. Therefore it is 
reasonable to ask for relations between them. We have chosen to show a connection between the 
multivariate umbrae they are related to, which can be translated in a connection between the 
t-th-order multivariate Bernoulli and Euler numbers. 

Lemma 3.9. If i is the multivariate Bernoulli umbra and 77 is the multivariate Euler umbra, 
then 2l = ^(77 — it) + 1. 

Proof. We have f(2v, z) = /(/,, 2z) = 2( Zl + ■■■+ z d )/[e 2 &+~ +z *> - 1] = f(n - u, §)/(*,, z). □ 

Theorem 3.10. (Relation between multivariate Bernoulli and Euler polynomials) We 

have 2^Bv^ (?r) = E[£^\x + t.L)] where l is the multivariate Bernoulli umbra. 

Proof. From Lemma 13.91 we have t.{2i) = t. [5(17 — u) + 1] = t.^(rj — u) + t.L, where last 

equivalence follows form the third equivalence in (|2.9p . The result follows since 2\ V \B$ (^) = 
E[{x + t.{2L)) v \. □ 
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Conclusions and open questions: multivariate time-space harmonic polynomials. 

A family of polynomials {P(x, t)}t>o is said to be time-space harmonic with respect to a stochastic 
process {X t } t > if E[P(X t ,t) | 3s] = P(X s ,s), for all s < t, where 3s = a [X T : r < s) is the 
natural filtration associated with {Xt}t>o- Recently [3] the authors have introduced a new family 
of polynomials which are time-space harmonic with respect to Levy processes and by which to 
express all other families of polynomials sharing the same properties. These polynomials are 
Appell polynomials and have the form E\(x + t.a) 1 ] for all positive integers i. By generalizing the 
definition of conditional evaluation given in [3] to the multivariate case, the multivariate Bernoulli 
and Euler polynomials should result to be time-space harmonic with respect to the multivariate 
Levy processes — t.i and ^[£.(w — I.T7)] respectively, whose probabilistic counterparts could be 
recovered via Propositions 12.8 1 and 12.91 Indeed, Corollary 13.51 shows that these polynomials share 
one of the main properties of time-space harmonic polynomials: when the vector of indeterminates 
is replaced by the corresponding Levy process, their overall mean is zero. We believe that the 
setting here introduced, together with the one given in [3], could be a fruitful way to build a 
theory of time-space harmonic polynomials with respect to multivariate Levy processes. 
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